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The superluminal propagation of neutrinos observed by OPERA collaboration can be interpreted 
as neutrinos traveling in a pseudoscalar potential which may be generated by a medium. The 
OPERA differential arrival time data set a constraint on the form of the pseudoscalar potential. 
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I. INTRODUCTION 

The OPERA collaboration recently reported the evi- 
dence of superluminal propagation of muon neutrinos [l| . 
An earlier report by the MINOS collaboration [H reached 
the similar conclusion with a lower statistical confidence. 
If the result is not due to an unknown instrumental effect, 
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this fact inevitably points towards new physics. Since the 
release of the OPERA data, many ideas have been dis- 
posed to interpret the data, e.g. The difficulties of 
various interpretations have been discussed One 
of the difficulties is to account for the energy dependence 
of neutrino speed anomaly A = {v — c)/c. While OPERA 
reports 
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A = (2.48 ± 0.28(stat.) ± 0.30(sys.) x 10"^ 
Ai = (2.18 ± 0.77(stat.) ± 0.30(sys.) x lO"'^ 
As = (2.76 ± 0.75(stat.) ± 0.30(sys.) x 10"''^ 



, E^^n GeV 
, E^^i = 13.9 GeV 
, E^2= 42.9 GeV 
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for the anomaly in the entire sample, Ei, < 20 GeV, and 
E^ > 20 GeV samples, respectively, SN 1987A neutrino 
arrival data poses a stringent constraint 

|A| < 2 x 10~^ E^ = 10 MeV . (4) 

The MINOS data, with less significance, suggests 

A = (5.1 ± 2.9) x 10-^ i?^ 3 GeV . (5) 

This further sharpens the A contrast in the 10 MeV to 
10 GeV range. 

Here we propose a superluminal neutrino theory with- 
out introducing Lorentz Invariance Violation (LIV) or 
particles with imaginary masses (tachyons) . We envisage 
a pseudoscalar potential under the influence of which the 
neutrino propagates from CERN to OPERA detector. 



II. SUPERLUMINAL NEUTRINOS IN A 
PSEUDOSCALAR POTENTIAL 

We introduce a pseudoscalar potential > 0, which 
can be energy-dependent, but is constant in space (or 
can be approximated as a constant for the spatial range 
in consideration). The Dirac equation for a neutrino with 
mass TOj/ in such a pseudoscalar potential can be written 
as 



[a.p + l3{m,y + 75^)] ip = E^ip, 



(6) 
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where a, /3 and 75 are the Dirac matrices. The natural 
units with c = 1 and h — 1 are adopted here. The 
neutrino wave function ip can be expressed in terms of 
two-component spinors u and v as 



i; = N 



(7) 



where N is the normalization constant. These two 
spinors u and v satisfy the coupled equations 



(a.p — (j))u — {E^ — m^)v ~ 
{a.p + (l))v ~ {El, + rn„)u — 



(8) 



Solving these two equations, we obtain the neutrino en- 
ergy 



For <C <IC E,^, one can write 

The group velocity of the neutrino is then 



Vn 



dE^ 
d\p\ 
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Ei. ' 
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where (j) — dcji/dEu, or 

A = w„ - 1 



2E?, E,, 



(9) 
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The condition of superluminal neutrinos (A > 0) can be 
written as 



III. 



CONSTRAINING THE FORM OF 
DATA 



(13) 
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The condition p3|) is satisfied for a wide range of 
forms. Below we attempt to construct an analytical form 
of (j) that gives superluminal motion and satisfies the ob- 
servational constraints. 
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C. Power law exponential 



The pseudoscalar potential of the form 



gives 



A 



a 

2 



Eo 



(18) 



(19) 



The superluminal condition is a < 1/2 + E^/Eq. How- 
ever, when Eu <^ Eq and E^, ^ Eq, the velocity differ- 
ence has the energy-dependence in the form of A oc E~^ 
and A oc E~^, respectively. Both dependences are too 
steep to account for the OPERA data in equations ([2]) 
and This form is therefore not favored. 



An energy-independent potential 
satisfies the condition p^ . so that 



0) naturally 



i.M 

E. 



2E?, 



> 1 
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Such a form is similar to many theories invoking LIV 
or tachyon particles e.g. [1, [l^j, which is ruled out 
by the OPERA differential speed data. The ratio of ve- 
locity difference in this model is A1/A2 — {£^^2/ E^^if . 
According to equations ([2]) and ([3]), the ratio between 
two average neutrino energies is about E^^2/E^^i 3.1. 
This corresponds to a predicted A1/A2 ~ 9.5, while the 
observed ratio is A1/A2 = 0.79lj!);5Q. So a more compli- 
cated form of (p with a much shallower E,^ dependence 
on A is needed. 



B. Power law 

For a pseudoscalar potential of the form 

^iE,)=AE:, (15) 

the superluminal condition (|13p can be translated to 

-0.5<a<0.5. (16) 

The ratio of velocity difference in this model is A1/A2 = 
{Eu^i/Ka?"'"^- Solving for a using the OPERA data, 
one gets 



0.72 < lal < 1.55 



(17) 



with the typical value \a\ 1.1. This range is incon- 
sistent with the condition (I16p . and is in the subluminal 
regime. Therefore the power law model cannot interpret 
the OPERA data. 



D. Power law logarithmic 

We consider the potential of the form 

'E: 



<I>{E,) 



\Eq J \Eq 
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The superluminal condition (|13l) can be translated to 



E, 

— >exp 
^0 



1 - 2a 



(21) 



for a < 1/2 (the branch a > 1/2 is not favored since the 
i?i/-dependence of A is very steep). The velocity differ- 
ence takes the form 
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E?, 
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(22) 

For a not very high Eo, the superluminal condition (j2ip 
is easily satisfied. The steep oc E~'^ dependence is com- 
pensated by other factors in equation so that a shal- 
low i?^-dependence can be achieved in the superluminal 
regime. By properly adjusting Eo and normalization C, 
the OPERA data can be interpreted. 



IV. DISCUSSION 

We have shown that if neutrinos travel in a pseu- 
doscalar potential, superluminal propagation is possible 
under the condition given in equation (jl3p . without vio- 
lating special relativity. In order to interpret the OPERA 
data, a shallow energy-dependence is required in the su- 
perluminal regime. A potential form similar to equation 
((20)) with a < 1/2 can meet such a requirement. 

The very smaU |A| derived for SN 1987A of 10 McV 
neutrinos is difficult to account for with such a potential. 
If one adjusts Eq to interpret OPERA data, the 10 MeV 
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neutrinos would be in the subluminal regime with a large 
|A| violating the data constraint. If one instead adjusts 
10 MeV to the transition point from super luminal to sub- 
luminal, i.e. El, E^_c — i^oexp( ^_^^^ ) ^ 10 MeV, then 
the energy-dependence in the 10 GeV range is too steep 
to satisfy the OPERA data. In order to incorporate the 
SN 1987A data, one needs to either argue for a local ef- 
fect of OPERA anomaly (e.g. the pseudoscalar potential 
is related to the density, gravity or magnetic fields in the 
neutrino propagation path in the earth crust) or admit 
that the potential (PO]) is not the correct form in the low 
energy regime. 

If, however, the pseudoscalar potential pO| can be ex- 
tended to the high energy regime, it is interesting to note 
that in the PeV energy range where the internal-shock- 
origin neutrinos from gamma-ray bursts (GRBs) are sup- 
posed to be generated [31 , neutrinos are still superlumi- 
nal but with a much smaller A ~ 10~^'^. For typical high- 



luminosity GRBs at cosmological distances, the lead time 
would be several years. For nearby low-luminosity GRBs 
[l^[lB], the lead time can be as short as several months. 
This can be in principle tested by independent discover- 
ies of a high-energy neutrino burst detected by IceCube 
and a later GRB detected in the same direction. This 
may be possible for a bright, nearb y G RB event such as 
GRB 030329. The TeV neutrinos [13, [Hi have a much 
larger A ^ 10~®, and hence, a leading time of decades to 
centuries, which is difficult to test observationally. 

Finally, we'd like to comment on that the normaliza- 
tion parameter C in ^U\\ adjusted to fit the OPERA data 
naturally leads to ^ mi/, but (j) ^ ttT's- So it is allowed 
that the Dirac equation with pseudoscalar potential ([6]) 
also applies to electrons and other spin 1/2 leptons, but 
those particles cannot be superluminal since their masses 
dominate the 4> term in the energy-momentum equation 

dSD- 
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